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Viscous fingers in a channel with surface tension anisotropy are numerically studied. Scaling 
relations between the tip velocity v, the tip radius p and the pressure gradient Px are investigated 
for two kinds of boundary conditions of pressure, when v is sufficiently large. The power-law 
relations for the anisotropic viscous fingers are compared with two-dimensional dendritic growth. 
The exponents of the power-law relations are theoretically evaluated. 

I. INTRODUCTION 

Growth patterns in Laplace fields and diffusion fields have been intensively studied. [H, Q In the problem of crystal 
growth from melt, dendrites grow in a diffusion field of temperature. A steadily growing dendrite is characterized by 
the radius of curvature p in the tip region and the growth velocity v. In the case of the absence of surface tension, 
only the product of p and v is determined as pv = 2Z)P(A) oc in two dimensions, where A is the supercooling and 
D is the diffusion constant. The surface tension anisotropy determines the value of p and v uniquely. The tip radius 
and the velocity satisfy the scaling relations v oc A^ and p oc A~^ in two dimensions. As a result, vp^ is constant, 
even when A is changed. [1, 01 

On the other hand, viscous fingers grow in the Hele-Shaw cell between two plates, where the pressure satisfies 
the Laplace equation. Q The pressure gradient is a control parameter in the problem of viscous fingering, which 
corresponds to the supercooling A in the problem of crystal growth. In a channel of finite width, Saffman and Taylor 
found a steadily growing solution of the form ^ 
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The effect of surface tension is not considered in the analysis performed by the Saffman and Taylor. In the case 
of the absence of surface tension, the finger width A cannot be determined uniquely. McLean and Saffman solved 
numerically steadily growing solutions of viscous fingers in the case of finite surface tension. Q The finger width is 
uniquely selected, even if anisotropy is absent. However, the finger width A approaches 0.5, when the tip velocity v 
and pressure gradient are increased to infinity. 

Dendritic patterns were observed even in experiments of viscous fingering in Hele-Shaw cells perturbed by etc hing 
one of their plates. [1, i, [M [H [H 

The etching induces some anisotropy. Matsushita and Yamada, [ll| and Couder |l2j | 
found the relation vp^ =const. for a single viscous finger in a certain parameter range. The scaling relation is the 
same as the dendritic growth in crystal growth. They considered that dendritic patterns in both crystal growth and 
viscous fingering have the same properties. Kessler et al. and Dorsey and Martin studied numerically anisotropic 
viscous fingering by using the solvability condition or solving a nonlinear eigenvalue problem. [Tsl. [T^ Kessler et al. 
predicted a scaling relation by WKB approximation as A ~ v^/'^ in the limit of large v. \'S\ The relations of curvature 
and velocity, or curvature and pressure gradient, were not reported in these simulations. 

In a previous study, we investigated some relations between p, v and the pressure gradient for a range of relatively 
small V [l^ using the original method developed by McLean and Saffman. Q In this study, we continue the numerical 
simulation and show scaling relations in the range of sufficiently large v, where some universal scaling relations are 
expected. The purpose of this study is to point out some difference between the dendritic growth in anisotropic 
viscous fingering and that in crystal growth from the scaling relations. We consider that the difference comes from 
the difference between the diffusion field and the Laplace field. 



II. MODEL EQUATION AND BOUNDARY CONDITIONS 

The model equation and boundary conditions are the same as those reported in the previous paper. We briefly 
explain them, and the details were written in refs.7 and 15. Viscous fingers with the constant pressure pq are pushed 
into a viscous fluid with the viscosity /x in a Hele-Shaw cell. The mean velocity u in the viscous fluid satisfies Darcy's 
law: u = — (6^/12/i)Vp = Vcj), where p is the pressure, is the velocity potential, and b is the spacing between the 
top and bottom plates of the Hele-Shaw cell. From the incompressibility of the fluid, satisfies the Laplace equation. 
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A viscous finger is assumed to grow in a channel of width 2a. Two kinds of boundary conditions for the pressure field 
p at the interface of the fingers have been studied. 

AP^^. (2) 

= ITT + + mfynimr"], (3) 

where T(9) denotes the anisotropic surface tension, w„ is the velocity normal to the surface, 9 denotes the angle of 
the tangent vector of the interface from the x-axis, and R is the radius of curvature on the surface, [p is the radius 
of curvature at the tip position.) The boundary condition (2) represents the Young-Laplace relation. The boundary 
condition (2) is equivalent to the Gibbs-Thomson effect in the problem of crystal growth, if the pressure is replaced 
with the temperature. Most authors, such as McLean and Saffman,[7] and Kessler et al.[l^ used this type of boundary 
condition for theoretical and numerical analyses. This boundary condition is useful for comparing the dendritic growth 
between the diffusion field and the Laplace field. The Young-Laplace relation is exactly satisfied, when the two fluids 
are stationary. Later, Park and Homsy derived the boundary condition (3) for the moving interface, when the ratio of 
the viscosities of two fluids is 0(Ca^/'^).[l^ Here, Ca is the capillary number Ca = fiv/T. This boundary condition 
is more realistic for the moving interface, however, rather complicated. Even this boundary condition is not available, 
when the capillary number (velocity) becomes too large. The boundary condition (2) is obtained, if the second term 
of eq. (3) is neglected, and (7r/4)r in eq. (3) is replaced with T. We study viscous fingers using the two boundary 
conditions and discuss the difference of the scaling relations caused by the boundary conditions. 

The dimensionless velocity vector {U, V) is expressed using the complex number a,s U — iV — ge^'^. The arclength 
S from the tip position along the interface is changed into another variable s satisfying < s < 1. At the finger tip, 
s = 1 and s = at a; = —oo on the interface. Because logg — i{6 — n) is an analytic function, Cauchy's integral 
theorem yields 

lo,im)--- f'-^ds'. (4) 
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As an anisotropy of the surface tension, the fourfold rotational symmetry of the form T{6) — To(l — e cos 46*) is 
assumed, where 9 — 6 — 71. The boundary condition (2) for the pressure is expressed as 

Kgs^ I (1 — e cos 4:9)qs = ^ cos 9, (5) 

where q — (1 — X)q, k = Tob'^ir'^ / {12 fiva^ {1 — \)^}. The boundary condition (3) is expressed as 
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where T{9) = (4/7r)To(l — ecos46') and k = Tob'^ir'^ / {12pLva?{l — A)^} are assumed for the coefficient of the first 
term in eqs. (5) and (6) to take the same form. The parameter g is defined as, g = 8a/(7r^6). The solutions to the 
coupled equations of eqs. (4) and (5) (or eqs. (4) and (6)) are numerically solved with the modified Newton method 
by discretizing the interface with N = 300 intervals. If q and 9 are obtained, the finger interface is expressed through 
the relation 

X + ly = / — as. (7) 
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III. SCALING RELATIONS FOR ANISOTROPIC VISCOUS FINGERS 



Figure 1(a) displays a viscous finger obtained numerically for e = 0.15 and 1/C — 12 fiv a^ / {Tob^ir'^) — 42843, which 
is a quantity proportional to the capillary number Ca = p-v /T and therefore the tip velocity v. The width of the finger 
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FIG. 1: (a) Viscous finger (solid curve) for e = 0.15 and 1/C = 42843, and the SafTman- Taylor solution (dashed curve) with 
A — 0.1115. (b) Finger width A as a function of 1/C for e — 0.05,0.1 and 0.15. (c) Double-logarithmic plot of 1/C vs A for 
e = 0.1 and 0.15. 



is A = 0.1115. The dashed curve denotes the Saffman- Taylor solution with A — 0.1115. The numerically obtained 
solution is well approximated using the Saffman- Taylor solution when A is small. Figure 1 (b) displays A as a function 
of 1/C for e = 0.05,0, 1 and 0.15. For e = 0, the finger width is always larger than 0.5. However, the finger width 
becomes smaller than 0.5 owing to the surface tension anisotropy. Figure 1(c) displays a double- logarithmic plot of 
1/C and A for e = 0.1 and 0.15 in the range of 1/C > 100. The dashed curve denotes a power of A = (l/C)"-^'*. A 
power-law relation appears in the range of large 1/C. The exponent 0.24 is very close to the exponent 0.25 predicted 
by Kessler et al. with WKB approximation. Figures 2(a)-2(c) display double- logarithmic plots of (a) 1/C vs the tip 
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FIG. 2: (a) 1/C vs p. (b) Pi vs 1/C. (c) Pi vs p for e = 0.4. The dashed lines in (a), (b), and (c) respectively denote power- law 
relations with exponents -0.49, 1.3, and -0.65. 

radius p, (b) Pi vs 1/C, and (c) Pi vs 1/C for e = 0.15, where Pi = A/C. The pressure gradient at a; ^ oo is a 
control parameter that corresponds to the supercooling A in the problem of crystal growth. The pressure gradient 
at X ^ oo is evaluated as —dp/dx = {12^/b'^)d4>/dx ~ {i2pLv\)/\P cx A/C. On the other hand, the tip velocity v is 
proportional to 1/C. Figure 2 therefore shows the power-law relations p cx l/v^ '^'^ , v oc P^'^, and p cx 1/P^-^^. The 
exponent of p vs 1/v is consistent with the scaling relation p oc 1/w" ''^ for the dendrite in the crystal growth. However, 
the other exponents are different from the scaling relation for the dendrite in the two-dimensional crystal growth: 
V cx and p oc 1/A^. The difference of the scaling relations comes from the difference between the diffusion field 
and the Laplace field. In the Laplace field, there is no characteristic length scale such as the diffusion length in the 
diffusion field. As a result, the pressure should decrease constantly with a constant pressure gradient Px in the range 
of X —>■ +O0 in an infinitely long channel. 

We can evaluate the exponents of the scaling relations. As shown in Fig. 1(a), the Saffman- Taylor solution is 
a very good approximation for the anisotropic viscous finger with a large tip velocity v. Then the tip radius p is 
proportional to (Px/dy'^ a.t y — for eq. (1), which is evaluated as A^/{7r^(l — A)}. For a sufficiently small X, p ^ A^. 
The prediction A cx (1/u)^''^ by Kessler et al. yields the relation p cx or vp^ =const. The pressure gradient 

—dp/dx is proportional to A/C cx Aw, which leads to the scaling relations: Pa: oc Aw oc v^^^ oc p~^l'^. That is. 
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and p oc Px These scaling relations are consistent with the numerical results shown in 
Figs. 1(c), and 2(a)-2(c) for e — 0.15. The scaling relations other than A oc (Xjv)^!^ have not yet been reported. The 
exponents evaluated in the previous study [l^ were different from the present values, because they were evaluated 
when the tip velocity v was not sufficiently large. 
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FIG. 3: (a) Viscous finger (solid curve) for e = 0.15, g = 10, and 1/C = 47156, and the Safi'man- Taylor solution (dashed curve) 
with A = 0.1445. (b) Double- logarithmic plot of the finger width A as a function of 1/C for p = 10 and 30 at e = 0.15. The 
dashed line denotes a power law with the exponent -0.18. 

Next, we investigate the case of the boundary condition (3). Figure 3(a) displays a viscous finger for e = 0.15, g = 10, 
and 1/C = 47156. It is also well approximated using the Saffman- Taylor solution. Figure 3(b) displays a double- 
logarithmic plot of A as a function of 1/C — \2pva? / [TqIP'tt'^) at e = 0.15 for g = 10 and 30. We find the power-law 
relation A oc 1/u" in the asymptotic region of 1/C — > oo. The exponent a is estimated at 0.18, which is smaller than 
the exponent a = 0.25 for the case of the boundary condition (2). Figures 4(a)-4(c) display double-logarithmic plots 
of (a) 1/C vs p, (b) Pi vs 1/C and (c) Pi vs p for e = 0.15 and g = 10. The dashed lines in Figs. 4(a)-4(c) denote the 
power-law relations p oc v oc P^''^^ , and p oc 1/P^ **, respectively. For the boundary condition (3), dominant 

terms are 

.(1 - A).9(-g.)3.8(l -ecosAff)^ / ^^.(-sing) ^ ^ _ ^^^^^ 

as \4Jo(1 — ecos40) / 

in the limit of large v. Since k oc 1/u, the main term on the left-hand side is evaluated as K^^^dd/ds. The order 
estimate by WKB approximation yields A^ ~ k^^^ or A oc 1/u^/^. Because the Saffman- Taylor approximation is good 
for a large v, the tip radius p oc A^ oc l/v^^"^. Then, the relation vp^ =const. is expected to be satisfied. The exponents 
for the scaling relations v oc Pi and p oc l/P^ are expected to be /3 = 6/5 and 7 = 2/5, respectively. The numerically 
obtained exponents a — 0.18, (3 = 1.21 and 7 = 0.44 are consistent with the order estimates a = 0.167, /3 = 1.2, and 
7 = 0.4. 



IV. SUMMARY AND DISCUSSION 



We have numerically studied steadily growing viscous fingers with anisotropic surface tension in the range of large 
1/C values. We have confirmed the relation vp'^ =const. even for anisotropic viscous fingering with the boundary 
condition (2); however, the exponents of the scaling relations v oc P^ and p oc l/PJ are respectively (3 — 4/3 and 
7 = 2/3 in the anisotropic viscous fingers, which are clearly different from the values /3 = 4 and 7 = 2 for the dendrites 
in two-dimensional crystal growth. Because the boundary condition is the same, the difference of the scaling relations 
comes from the difference between the diffusion field and the Laplace field. 

If the boundary condition (3) is assumed, the relation vp'^ =const. is not satisfied even in the limit of large v; 
however, vp^ =const. with 5 ^ 1/(2q;) ~ 3 is expected to be satisfied, because the velocity effect dominates in the 
boundary condition in the limit of large v. 

In the experiments of Matsushita and Yamada, and Couder, the relation of vp^ appears in a relatively small range 
of V. They have also observed that the tip radius p becomes nearly constant in the range of large v. Matsushita and 
Yamada suggested that the crossover originates from the second term including the v of the boundary condition (3). 
We have found the relation vp^ =const., but our numerical simulation is not always consistent with these experiments. 
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FIG. 4: Double-logarithmic plots of (a) 1/C vs p for p = 10 and e = 0.15, (b) Pi vs 1/C for g = 10, and (c) Pi vs p. The 
dashed lines in (a), (b), and (c) respectively denote power-law relations with exponents -0.36, 1.21, and -0.44. 



in that we have not found p =const. in the range of large v for the boundary condition (3). In the experiment of 
Matsushita and Yamada, a dendrite grows along a straight groove. In the experiment of Couder, a bubble is attached 
to the tip region of a viscous finger. We consider that the boundary conditions and experimental setting of these 
experiments are different from those in our numerical simulations. 

In the experiments by Ben- Jacob et al. Q, McCloud and Maher and Honda et al. [l3|, many straight grooves 
are etched on the plates of the Hele-Shaw cell. Our numerical simulation corresponds to these experiments when v 
is sufficiently small, where the surface tension dendrites appear. However, in these experiments, the surface tension 
dendrites disappear, instead, kinetic dendrites that grow along the direction of the groove appear in the range of large 
V. It is because the spacing b in the groove regions is larger than that in the other regions, and therefore the mean 
velocity is larger owing to Darcy's law: u — ~{b^ /l2fi)Vp. Some kinetic effect of the large value of v and the effect 
of the uniform spacing b are considered in the boundary condition (3). However, the effect of the heterogeneity of 
the spacing b or the essential three-dimensional effect should be considered to discuss the kinetic dendrites growing 
along the direction of the grooves. We are not sure whether the effect of the heterogeneity can be considered by the 
modification of the boundary condition through a 0-dependent parameter, such as b{9), or direct three-dimensional 
simulations are necessary. 

To summarize, we have numerically shown that the scaling relations change with the fields, in which patterns grow, 
and the boundary conditions. Our results are mathematical ones, and the relevance of the scaling relations to actual 
experiments remains insufficient now as discussed above. We hope that the scaling relations will be confirmed by 
suitable experiments 
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